Study of the optimal conditions for NV center formation in type lb diamond, using 
photoluminescence and positron annihilation spectroscopies. 
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We studied the parameters to optimize the production of negatively-charged nitrogen-vacancy 
color centers (NV~) in type lb single crystal diamond using proton irradiation followed by thermal 
annealing under vacuum. Several samples were treated under different irradiation and annealing 
conditions and characterized by slow positron beam Doppler-broadening and photoluminescence 
(PL) spectroscopies. At high proton fluences another complex vacancy defect appears limiting the 
formation of NV~. Concentrations as high as 2.3 x 10^* cm~"^ of NV~ have been estimated from PL 
measurements. Furthermore, we inferred the trapping coefficient of positrons by NV~. This study 
brings insight into the production of a high concentration of NV~ in diamond, which is of utmost 
importance in ultra-sensitive magnetometry and quantum hybrid systems applications. 

PACS numbers: 81.05.ug, 61.72.jn, 61.72.jd, 61.72.Cc, 78.70.Bj 



I. INTRODUCTION 



II. EXPERIMENTAL PROCEDURE 



Owing to its unique properties of perfect photolu- 
minescence stability and optically detectable electronic 
spin resonance with a long coherence time, the nitrogen- 
vacancy (NV~ ) color center in diamond has been the ob- 
ject of a large number of studies over the past few yearsi*^ 
from quantum information processing, to highly sensitive 
magnetometry'^'^ and cellular imaging^. Most of these 
studies rely on the properties of isolated single centers, 
but some applications can benefit from a high density 
of NV~ centers while keeping a long spin coherence — , 
e.g. the realization of collective quantum memories^, hy- 
brid quantum circuits in which superconducting qubits 
are coupled to NV~ electron spins^, as well as ultrasensi- 
tive magnetometry at the micrometer scale^*^. Further- 
more, the application of diamond nanoparticles as labels 
for bio-imaging also requires a high NV~ contenfeiSiii. 
A well-known route to produce a high concentration of 
NV~ centers is the irradiation by high-energy particles 
and the subsequent annealing of nitrogen-rich diamond. 

In this work we determined the optimal parame- 
ters for the production of such a high NV^ concen- 
tration. Several type lb diamond samples were irra- 
diated with 2.4 MeV protons at fluences ranging from 
10"'^^ to 10^^ cm^^, and annealed under vacuum at tem- 
peratures from 600°C to 1000°C for durations from 1 
to 20 hours. Each sample was characterized at differ- 
ent stages of the NV^ formation process (before/after 
irradiation and after annealing) using photolumines- 
cence (PL) spectroscopy and slow positron beam based 
Doppler-broadening spectrometry, referred in short as 
Positron annihilation spectroscopy (PAS). The latter 
technique provides information about the nature of the 
vacancy defects produced along the depth probed by the 
positronsi^ii^. 



The samples used are one-side polished 3x3x0.5 mm"^ 
type lb HPHT diamond plates (Element Six Ltd.) with 
a nitrogen concentration specified below 200 ppm (3.5 x 
10^^ cm"'^). Sample irradiation was performed with a 
2.4 MeV proton scanning beam from a Van de Graff ac- 
celerator. During irradiation, the temperature of samples 
was maintained below 80° C using a water-cooled sample 
holder. SRIM simulationsAi show that such protons cre- 
ate an almost constant concentration of atomic displace- 
ments in the first 25 fim beneath the surface and stop at 
about 35 /im producing cascades, which increase sharply 
the local damage concentration. The sample annealing 
was performed under vacuum. PL measurements were 
done with a confocal microscope equipped with an imag- 
ing spectrograph. We estimate the NV~ concentration 
by quantitative comparison of the emission spectrum to 
the one of a single color center 



For PAS measurements we used an accelerator deliv- 
ering monoenergetic positrons with energy varying from 
0.1 to 25 keV (for a complete description see Ref.^'*). A 
Doppler broadening spectrometer coupled to the positron 
accelerator allows us to measure the energy of the gamma 
photon emitted from the electron-positron annihilation 
and hence to probe the vacancy defects. We determined 
the usual PAS parameters: S (resp. W) corresponds to 
the annihilation fraction of low (respectively high) mo- 
mentum electron-positron pairs, probing predominantly 
valence (respectively core) electronsi^. S and W change 
with an increase in the vacancy defect concentration or 
with a change in the nature of the defects probed. The 
measurements of S{E) and W{E) as a function of the 
positron energy E allow to probe the sample as a func- 
tion of the depth in the first 5 /im under the surface. 
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III. RESULTS AND DISCUSSIONS 



PAS and PL measurements were first performed on un- 
treated as-received diamond samples, then on irradiated 
but not annealed samples, and finally on irradiated and 
annealed ones. 

For all PAS measurements presented in this work, the 
positron annihilation characteristics S{E) and W{E) as 
a function of the positron incident energy can be 
viewed as resulting from a combination of surface and 
bulk characteristics. S{E) and W{E) exhibit plateaus 
at high positron energy (12-25 keV) (see Fig. [irb)-(c) 
and inset of Fig. HJa)), indicating that the contribution 
of the surface can be considered as negligible at the cor- 
responding depths > 1 and hence {S, W) values are 
solely those of the bulk. From now on, we shall only refer 
to these plateau values and we consider that the distri- 
bution of vacancies is homogeneous in the region probed 
by the positrons. 



A. Photoluminescence and positron annihilation 
measurements done on as-received samples 



Figure [TJa) shows that the as- received sample exhibits 
a sharp PL peak at around 573 nm corresponding to the 
first order Raman peak of pure diamond at the excitation 
wavelength of 532 nm. A broader peak around 690 nm 
indicates the presence of native NV~ color centers. A 
concentration of 8.3 x lO^'^ cm"'^ of native NV~ is esti- 
mated from the integration of this peaki^. 

PAS measurements conducted on four different un- 
treated samples show similar S{E) and W{E) graphs 
(Fig. [IKb)-(c)), indicating the homogeneity of diamond 
quality in the set of samples with respect to the atomic 
structures to which the positrons are sensitive. From 
these measurements, we inferred: 5ref = 0.351 ± 0.001 
and Wrcf = 0.0834 ±0.0006. 

A fit of the data using VEPFIT^^ allows us to deduce 
a positron effective diffusion length, of about 150 nm. 
This value is comparable to the value of 111 nm, which 
was found for in type lb diamond^^. This shows 
that we can consider the bulk annihilation characteris- 
tics (S'rofjWrcf) measured in the as-received samples as 
those of the pure diamond lattice (S'l, VFl)- This implies 
that the concentration of vacancy-related defects present 
in the untreated type lb samples is below the sensitivity 
threshold of the PAS technique corresponding to a value 
of about 10^5 cm-3 [Fig. (IKb)]. 

Moreover, in the following, we will only consider the 
relative values S'/S'l and W/Wi, which are less sensitive 
to the experimental setup. 
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FIG. 1. (a) PL of as-received diamond sample (cw laser exci- 
tation at the wavelength of 532 nm, laser power: 13.5 /^W), 
(b) Low momentum annihilation fraction S and high momen- 
tum annihilation fraction W of positrons as a function of their 
incident energy E, in four as-received diamond samples iden- 
tified by references DlbES022-25. 



B. Photoluminescence and positron annihilation 
measurements done after irradiation and prior to 
thermal annealing 

1. Photoluminescence measurements and inferred NV~ 
color centers concentration 

PL measurements done on samples irradiated at dif- 
ferent proton fluences [Fig. [2jb)] display a broadband 
around 690 nm due to NV~ photoluminescence and a 
higher energy PL band cut below 550 nm by the dichroic 
beamsplitter. This band is barely visible at the low- 
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FIG. 2. PAS and PL measurements realized on diamond sam- 
ples after proton irradiation (fluences from lO^'^ to 10^^ cm~^) 
and before annealing. {a){S/ S-l,W/W-l) scattered plot, with 
S and W taken at their plateau values in the 12-25 keV en- 
ergy range (min of 10 measurements per fluence), Sl and Wl 
being the diamond lattice parameters'^; inset: S and as a 
function of the positron energy (0.1-25 keV). (b) PL spectra 
measured on the same samples. Excitation with a cw laser 
(wavelength: 532 nm, power: 13.5 /xW). The narrow peak at 
573 nm is the first order Raman scattering. 



est fluence but becomes predominant at the highest flu- 
ences. It is tentatively attributed to the 113 (N-V-N) cen- 
ter emission^^. The absence of the GRl center emission 
(around 741 nm) associated to the neutral monovacancy 
is noteworthy. It has been shown that in type lb dia- 
mond the monovacancy is preferentially in the negative 
charge state^i due to the presence of the nitrogen impu- 
rity, which can donate an electron to the monovacancy. 
Both charge states coexist when there are not enough ni- 
trogen donors with respect to vacancies. For the fluence 
5 X 10^^ cm^^ the maximal vacancy concentration can 
be estimated by that of the atomic displacements cal- 
culated by SRIM^- to be 5 X 10^^ cm""^, which is larger 
than the specifled nitrogen concentration 3.5 x 10^^ cm~'^. 
The absence of the GRl line therefore indicates that the 
real monovacancy concentration is lower than the one of 
atomic displacements, as a result of vacancy-interstitial 
recombination or non-photoluminescent complex vacancy 
formation. 
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FIG. 3. Concentration of NV~ in different diamond samples 
irradiated at 5 x lO'^ cm~^ as a function of the annealing 
temperature (600-1000°C) and duration (1 or 10 hours). 



2. Positron annihilation measurements 



Figure [2] (a) shows that when we increase the fluence, 
5/S'l increases and W/Wi, decreases indicating the for- 
mation of vacancy defects. For fluences < 5 x 10^"^ cm^^, 
the data are well fltted by a straight line Dl starting 
from the as-received sample measurements. This line is 
the signature of an annihilation state additional to the 
lattice one and consecutive to irradiation. Since the an- 
nihilation fraction in this state increases with the fluence, 
we suggest that it corresponds to the negatively charged 
monovacancy, which is the most simple defect that can 
be created in this type of nitrogen-rich diamondiS. 

For fluences > 5 x 10^^ cm~^ the {S, W) points are 
located above Dl, indicating that a new type of vacancy- 
related defect appearsi^. This claim is supported by the 
fact that we measure larger relative changes of S val- 
ues than the one of the neutral or negative monovacancy 
in diamond^-^ (5'v°/S'l = 1.08 ± 0.01 and S'v-/5'l = 
1.079 ± 0.002) indicated by the horizontal black line in 
Fig. [5] (a)^ from the intersection of this monovacancy 
horizontal line with Dl. In order to check if NV~ cor- 
responds to this defect center, we estimated its highest 
concentration from the PL spectrum [Fig. [Hb)], to be 
1.5 X 10^^ cm-3 at fluence 5 x lO^'^ cm"^. Such a NV" 
concentration is hardly detectable by PAS as proved by 
further measurements on annealed samples [Fig. [Sjb)]. 
Moreover this concentration is over-estimated due to the 
contribution of the high-energy band in the PL spectrum. 
Therefore NV~ centers present in the samples after irra- 
diation but prior to annealing can only marginally con- 
tribute to the (S'/S'l, VF/Wl) shift. This deviation to the 
Dl observed at the highest fluences is most likely related 
to multiple vacancy defects. 
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FIG. 4. PL spectra (a) and inferred NV~ concentration (b) 
in diamond samples as a function of the proton fluence, under 
identical 800° C and 20 h annealing conditions. At high pro- 
ton fluence, nitrogen converts also into the non-paramagnetic 
neutral NV° cente^^l characterized by a zero-phonon line 
(marked by **) at 575 nm wavelength, while the one char- 
acteristic of the NV~ is at 637 nm (*). In graph (b) the 
dashed (dotted) line displays the concentration of atomic dis- 
placements predicted by SRIM (resp. [N]=200 ppm). 
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FIG. 5. PAS measurements on samples irradiated at different 
fluences and annealed at 800°C for 20 h. (a) {S/Sl,W/Wl) 
scattered plots. Filled (open) symbols refer to irradiated but 
non-annealed (resp. annealed) samples. Black squares refer 
to the untreated sample, (b) S/Sl vs. NV~ concentration in 
annealed samples; plain curve: best fit by relation 



NV center formation for fixed annealing conditions 
( 80(f C, 20 hours ) and various fluences 



C. Post-thermal annealing studies 

1. NV~ center formation at fixed proton fluence of 
5 X lO^'' cm"^ and various annealing conditions 

To investigate the effect of annealing time and tem- 
perature, diamond samples were irradiated at a proton 
fluence of 5 x 10^^ cm~^ and annealed at temperatures 
ranging from 600°C to 1000°C, during 1 or 10 hours. 
Figure [3] displays the estimated NV^ concentration for 
the different conditions. A significant increase in the 
NV~ concentration above 700°C is observed. Further- 
more [NV^] saturates at ~ 1.7 x 10^^ cm"^. Note that 
from 800°C to 1000°C, the difl?erence in NV" concentra- 
tion for the two annealing times becomes smaller, even- 
tually vanishing at 1000°C. The NV~ formation kinetics 
is thus enhanced upon temperature increase. 



We also studied the influence of the irradiation fluences 
(from 10^2 cm-2 to lO^"^ cm ^ ) on NV center formation 
for a given annealing condition (800°C, 20 hours). The 
PL spectra of FigUfa) show an intensity increase with 
the fluence up to 5 x 10^^ cm"^. Yet at all the fluences, 
the inferred NV~ concentration plotted on Fig. Hfb) is 
one order of magnitude lower than that of the atomic dis- 
placements predicted by SRIM, probably due to vacancy- 
interstitial recombination. The maximum NV~ concen- 
tration reached is 2.3x 10^® cm~^. Note that a lower NV"~ 
concentration is measured for the highest fluence. This 
decrease could have been attributed to the onset of di- 
amond amorphization but confocal Raman spectroscopy 
did not reveal any sp2 carbon signaliS,. However, PAS 
measurements provides an explanation as shown in the 
next paragraph. 

Figure ^a) displays S{W) for the different fluences 
before and after annealing. At all the fluences, we ob- 
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serve a decrease of S concomitant with an increase of 
W after annealing, indicating that the vacancy defect 
concentration diminishes and/or the nature of the de- 
fects changes. {S/ Sh,W/Wh) points after anneahng are 
aligned on a new straight line D2 except the highest flu- 
ence one. Since PL measurements evidence the creation 
of NV~ defect, we propose that it is the defect associ- 
ated to D2. This claim is supported by previous work 
showing an increase of W due to nitrogen decoration of 
vacancies in diamond32,. Note that for 5 x lO^'* cm~^ and 
5 X 10^^ cm~^ fluences {S/Si,,W/Wi^) points which are 
initially off the monovacancy Dl line finally end onto D2 
after annealing. It could be due to the positron trapping 
rate balance between different vacancy traps shifting the 
S and W values to the annihilation characteristics of the 
NV~ vacancy defect which has a high trapping coeffi- 
cient, without excluding that annealing can partially re- 
move the complex vacancies. In contrast, for the highest 
fluence the {S/Sl, W/Wl) points almost did not move af- 
ter annealing, and still fall off the D2 line, which indicates 
that there is a significant fraction of complex vacancies 
that do not evolve into NV~ centers, thus limiting the 
NV~ creation yield. It would be worthwhile to anneal at 
temperatures higher than 1000°C. This strategy was suc- 
cessful in increasing the NV~ electronic spin coherence 
as a result of paramagnetic defect destruction^^ including 
divacancies associated ones like R4 or W62i. 



3. Trapping coefficient of positrons by NV suggested by 
positron annihilation and photolummescence measurements 

The combination of PAS and PL measurements also 
allows us to estimate the trapping coefficient of positrons 
by NV~ . Fig. [5ljb) displays S* as a function of the NV^ 
concentration inferred from PL spectra [Fig. |4l[a)]. The 
data can be fitted by a relationship derived from rate 
equations considering the lattice and vacancy defects 
(mostly NV^) as the positron trapping sites: 

S Al ^ A^NV- C/ Cl, 5'nv- ^^-^ 

Sl Al + /iNv-C'/CL Al + /iNV- C'/Cl S'l 

where Al is the lattice annihilation rate, C and Cl = 
1.8 X 10^'^ cm~'^ the NV~ center concentration and di- 
amond lattice atomic density respectively, /^nv- is the 
associated positron trapping rate. Taking Al — 10 ns^^ 
from Refii^, the best fit to the data is obtained for 
Sj^y-/Sl = L074 ± 0.005 and /xnv- = (7.6 ± L6) x 
10^^ s~^. The same type of fit can be performed for 
W/Wl, yielding Wj^y- /Wl = 0.82 ± 0.05 and /.tNv- = 
(6.7 ± 1.8) x lO^s s-i. The two values of /^nv- ob- 



tained are consistent, and we take their mean value 
/ijvjy- = (7.2 ± 1.8) X 10^^ as the one inferred from 
PAS. In silicon the positron trapping coefficient has been 
found to be (0.8 ± 0.4) x 10^^ s'^ for the neutral di- 
vacancy25 and 1.8 x 10^^ s^^ for the (VP)" defect^l. 
The more efficient trapping of V" in diamond compared 
to (VP)~ in silicon can be explained by: (1) a smaller 



TABLE I. Positron annihilation parameters and trapping rate 
fi for the vacancy defects in diamond, as estimated in this 
work (except for Sy-fSi,). 



Positron parameters V 


NV" 


S/^L (Sl = 0.351 ± 0.001) 1.079 ± 0.003^^ 


1.074 ±0.005 


W/Wi,{Wl = 0.0834 ± 0.0006) 0.82 ± 0.05 


0.87 ±0.01 


(xlO^'' s-i) 


7.2 ± 1.8 



static dielectric constant (5.6 for diamond and 11.6 for 
silicon), hence a longer range coulomb potential in the 
case of diamond, (2) a binding energy of 0.3 eV for the 
Rydberg state of the V~ electron, larger than the one of 
0.1 eV in silicon, leading to a significant reduction of the 
positron detrapping in diamond. Table |T] summarizes the 
PAS data related to the defects detected in this work. 



IV. CONCLUSION 

We have optimized the formation of a high-density of 
NV~ centers in type lb diamond by proton irradiation 
and subsequent thermal annealing. We showed that the 
conversion efficiency of nitrogen into NV~ increases with 
the irradiation fiuence, up to a point where further in- 
crease leads to a reduction in this conversion rate, most 
likely due to the creation of complex vacancies stable 
under 800° /20 h annealing. We achieved the maximum 
NV~ concentration of 2.3 x 10^® cm~'^ (13.5 ppm), cor- 
responding to a conversion efficiency [NV~]/[N] of about 
6.7% ([N]~ 200 ppm). The positron trapping coefficient 
by NV~ centers in diamond was estimated for the first 
time. 
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